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Abstract 

A previous model where the low-lying axial-vector mesons are dynamically gener- 
ated, implementing unitarity in coupled channels in the vector-pseudoscalar (VP) 
meson interaction, is applied to evaluate the decay widths of the ai(1260) and 
6i(1235) axial-vector mesons into irj. Unlike the case of the a\, the h\ radiative 
decay is systematically underestimated at tree level. In this work we evaluate for 
the first time the loop contribution coming from an initial VP vertex. Despite the 
large superficial divergence of the loops, the convergence of the relevant loops can 
be established by using arguments of gauge invariance. The partial decay widths 
obtained agree with the experimental values within uncertainties, and we show that 
the loop contribution is crucial in the b\ case and also important for the a\ case. 



1 Introduction 

The unitary extensions of chiral perturbation theory (%PT) have allowed to extend the 
range of energies where the hadron interaction can be studied. At the same time they have 
also shown that many meson and baryon resonances are dynamically generated and can 
be interpreted as quasibound states of pairs of hadrons in coupled channels [1]. A case 
very well studied is the one of the interaction of the octet of pseudoscalar mesons [2-6] 
from where a nonet of scalar mesons are generated. Much less studied is the case of the 
interaction of vector mesons with pseudoscalar mesons, where two independent works [7, 8] 
have shown that the axial vector mesons can be generated dynamically. This novel idea 
should be confronted with experiment to test the accuracy of its predictions. Some of these 
predictions have already been tested in Ref. [9]. Contrary to other pictures like quark 
models, where external sources are coupled to the quarks, in the dynamically generated 
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picture one assumes that the largest weight of the wave function is due to the two meson 
cloud, and consequently, the coupling of external sources proceeds via the coupling to the 
meson components. One interesting test which brings light into this issue is the radiative 
decay of the resonances. This is the purpose of the present work where we concentrate 
on the radiative decay of the bf and af axial vectors into The af radiative decay 

has been studied within different contexts, for instance vector meson dominance is used 
in [10,11], relating the radiative decay with the /wr decay of the af. Chiral Lagrangians 
with vector meson dominance (VMD) are also used in [12] to obtain the radiative width of 
af — > 7r + 7. An SU(3) symmetric Lagrangian is used in [13] to account for strong decays 
of the axial vector mesons and by means of VMD the amplitude for af — > 7r + 7 is studied 
and related to the one of [12]. A common feature of these works is that the bf — > 7r + 7 
reaction is not discussed and its evaluation in [13] using VMD along the same lines as the 
af — > 7r + 7 gives rise to a decay rate substantially smaller than experiment. The bf — > 7r + 7 
decay is also neglected in the analysis of [11] citing the small rates obtained. 

The rates of af — > 7r + 7 and bf — > 7r + 7 are also evaluated in [14] using quark models for 
the ai — > irp and bi — > ixoj and VMD to relate these amplitudes with the radiative decay. 
It is emphasized there that because of the factor 1/3 of the c^7 coupling relative to the 
one of p7 there is a reduction factor of 1/9 for the radiative decay bf — > 7r + 7 compared to 
that of the af — > 7r + 7 decay, resulting in a ratio of these two rates in contradiction with 
experiment (this is the same argument found in [13] as responsible for the small rate of the 
bf — > 7r + 7 decay). 

In the present work we shall also use the tree level VMD amplitudes, but in addition, 
the nature of the axial vector mesons as dynamically generated resonance provides a strong 
coupling to K*K and K*K, and subsequent loops with these intermediate states and the 
photon emitted from these constituents should be considered. We show that the loops are 
very important, particularly for the case of the bf — > 7r + 7 decay, and the simultaneous 
consideration of the VMD amplitudes at tree level and loop contributions leads to a good 
description of both radiative decays. 

We shall also show some technical details involving loops with vector mesons. Using 
arguments of gauge invariance and the Feynman parametrization, one can prove that the 
loops involving one vector meson and two pseudoscalar mesons are finite, in spite of the 
large degree of superficial divergence. This was found in [15-18] with the loops involved in 
radiative decay of the <p containing three pseudoscalar mesons. 

2 Formalism 

In Ref. [8] most of the low-lying axial vector mesons were dynamically generated from the 
s-wave interaction of the octet of vector-mesons with the octet of pseudoscalar-mesons by 
using the techniques of the chiral unitary theory. With the only input of a chiral Lagrangian 
for a vector and pseudoscalar (VP) mesons and the implementation of unitarity in coupled 
channels, these resonances show up as poles in the second Riemann sheet of the unitarized 
scattering amplitudes. By evaluating the residues of the scattering amplitudes at the 
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pole positions, the couplings of the dynamically generated axial-vector resonances to the 
different VP channels can be obtained. By using these couplings, we found a nice agreement 
with the experimental VP partial decay rates, despite the fact that no parameters were 
fitted to experimental data of the axial-vector mesons. 

In view of the dominant contribution of the VP channels in the building up and decay of 
the axial-vector resonances, our starting point to study the radiative decay of the b\ and a\ 
is the transition of these resonances into the allowed VP channels and attaching the photon 
to the relevant meson lines and vertices. The first kind of mechanisms considered are the 
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Figure 1: Feynman diagrams contributing to the radiative axial- vector meson decay. 



tree level vector meson dominance (VMD) contributions, shown in Fig. [Tb). Furthermore, 
the radiative decay can also proceed from loops of the VP pair with the photon emitted 
from either the pseudoscalar or the vector meson leg, Fig. CD 1 )) and c). A diagram with 
the photon directly emitted from the VPP vertex is needed to ensure gauge invariance, 
but we will explain later on (after Eq. ( BSD ) that, using arguments of gauge invariance, we 
do not need to evaluate it directly. On the other hand, another kind of loops containing 
the VPj and VVP vertices is also possible, however, they involve two abnormal intrinsic 
parity vertices, and hence its contribution should be rather small compared to those already 
considered. This is indeed the case in the analogous loops present in the radiative decay 
of the meson, as it was found in [19]. 

The intermediate VP states possible in the loops are those used in Ref. [8] to build up 
the axial vector mesons. These are, for the b\: l/y2(K*K + K*K), 07r, ujtc, pr), and for 
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the a\.\j\/2[K*K — K*K), pir. Note however, that the coupling of the 07r, luh, pr\ to the 
final pion violates G— parity and hence these channels do not contribute to the b\ radiative 
decay. Thus, only the diagrams in Fig. [T] must be evaluated. 

Let us start with the evaluation of the tree level contributions. For the V7 vertex we 
use the amplitude 

£y 7 = -e\ v F v M v e v ■ e (1) 

with Ay = 1, 1/3, -V2/3 for p, u and respectively, F v = 156 ± 5MeV [19], My is the 
vector meson mass and ey and e are the vector-meson and photon polarization vectors 
respectively, and e is taken positive. 

The axial-vector meson coupling to VP can be expressed [8] as 

tAVP = 9AVP £a ■ ey (2) 

where €a is the axial- vector meson polarization vector. The couplings gAVP are obtained in 
Ref. [8] by evaluating the residues at the poles of the VP unitarized scattering amplitudes 
and are given in table VII of that reference. Note that in Ref. [8] the couplings are given 
in isospin base and for given G— parity states, hence the appropriate projection to charge 
base has to be done. In Ref. [8] no theoretical errors were quoted for these AVP couplings. 
However, for the purpose of evaluating the theoretical uncertainty in the calculations of 
the present work, we have estimated the uncertainties in these couplings in the following 
way: for the b\ case, we have considered the change in the couplings due to a resonable 
uncertainty of 10% in the only free parameter of the model, the subtraction constant 
a ~ —1.85 (see Ref. [8] for details). We consider further uncertainties from changing /, 
as will be explained after Eq. ([7j). For the a\ case, in Ref. [8] the mass obtained was 
1011 MeV, somewhat below the nominal mass in the PDG [20], 1230 MeV. (Note however 
that the total width is 250 — 600 MeV in the PDG, which gives an idea of the uncertainty 
in the mass). In Ref. [8], the mass was obtained with a value of a = —1.85. If we use the 
value a = —1.1 and / = fx, we obtain a mass closer to the nominal of 1080 MeV, and it 
is not easy to get larger mass. In this case the coupling to pit, the dominant channel, is 
increased by ~ 25%. From there, we get an idea of the uncertainties in the a% couplings. 

This leads to the following radiative decay amplitudes for the tree level diagrams, 
(Fig. Eh)): 




(3) 

In the evaluation of the loops an apparent problem arises given the large superficial 
divergence due to the loop momentum dependence of the vertices and the q^q u /My terms 
of the vector meson propagators. However, we will explain in detail how one can circumvent 
this problem invoking gauge invariance and using a suitable Feynman parametrization of 
the loop integrals. 
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Since the only external momenta available are P (the axial-vector meson momentum) 
and k (the photon momentum), the general expression of the amplitude can be written as 



T = (4) 

with 

= ag^ + b P»P V + c P*k v + d WP V + e k"k u (5) 

Note that, due to the Lorentz condition, e^P^ = 0, e v k v = 0, all the terms in Eq. (j4]) 
vanish except for the a and d terms. On the other hand, gauge invariance implies that 
T^ v k v = 0, from where one gets 

a = -dP-k. (6) 

This is obviously valid in any reference frame, however, in the axial-vector meson rest 
frame and taking the Coulomb gauge for the photon, only the a term survives in Eq. (jl]) 
since P = and e = 0. This means that, in the end, we will only need the a coefficient 
for the evaluation of the process. However, the a coefficient can be evaluated from the 
d term thanks to Eq. ([6]). The advantage of doing this is that there are few mechanisms 
contributing to the d term and by dimensional reasons the number of powers of the loop 
momentum in the numerator will be reduced, as will be clearly manifest from the discussion 
below. 

Let us start by evaluating the diagram []J) with the photon emitted from the pseu- 
doscalar leg for the bi — > 7r + 7 decay channel (the other channels are totally analogous). 
We will call this diagram type-b, in contrast with the type-c, with the photon attached to 
the vector- meson leg which will be evaluated later on, (see Fig. [1]). 

For the evaluation of this diagram we also need the VPP and PP7 vertices. The VPP 
Lagrangian used is (see Ref. [21] for normalizations used) 

£ VPP = -iJL(V>[d ll P i P]), (7) 

where and P are the usual SU(3) matrices containing the vector and pseudoscalar 
mesons. In Eq. flTJ) (...) means SU(3) trace and g = —MyGy/f 2 , where My is the vector 
meson mass, Gy = 55 ± 5 MeV [19] and / is the pion decay constant that we take from 
93 MeV to 1.15 x 93 MeV to take into account the uncertainty due to the use of / instead 
of fx which could actually enter some of the expressions. These uncertainties in the 
parameters, together with the uncertainties in other couplings of the theory, will be taken 
into account later on in the evaluation of the theoretical uncertainties of our results. 
The PP7 vertex is easily obtained from the lowest order ChPT Lagrangian [22] 

£ = ^(D^D»U) (8) 

where the photon field appears in the covariant derivative. 

The amplitude for the type-b mechanism (see Fig. [Tj?) for the bf — > 7r + 7 reads: 
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with 



d 4 q 



(27r) 4 q 2 — m\ (q — k) 2 — m\ (P — q) 2 — m 2 K * 
[ (P _ q) . (2k - 2P) + {P- q) 2 ] 



-(2k - g) M - <v 



2q u . 



(10) 



In Eq. ([9]), gb 1 K*K is the coupling of the bi to the K*K and K*K G-parity positive com- 
bination as defined in Ref. [8]. By looking at Eq. (flOl) one can see that in the worst case 
the loop integral, as it is written, is quadratically divergent. At this point, we can take 
advantage of the fact that we only need to evaluate the contribution to the d term of 
Eq. ([5]), as explained above. In fact, the most divergent term, the one with (P — q) 2 , does 
not contribute to the that term. Indeed, we can write 



(P - q) 2 = (P- qf - m 2 K , + mj 
and then the two first addends of the right hand side of the above equation give 



11^ 



d 4 q (P - qf 



m 



K" 



1 



(27r) 4 (P — q) 2 — m 2 K » q 2 — m 2 K (q — k) 2 — m 2 K 



q^v 



12) 



which does not depend explicitly on P, since the (P — q) 2 — m 2 K * cancels the propagator 
where the P appears and, therefore, this integral cannot give contribution to the d term. 
Hence, for the purpose of evaluating the fc M P, contribution, Eq. fflOl can be written as 



d\ 



(2vr) 4 q 2 - m\ (q - k) 2 -m\(P- qf 



m 2 K , 



-4fc„ + 2q p 



1 - 



(P - q) ■ (2k - 2P) 



m 



K* 



ni 



K* 



q v 



(13) 



which has one power less in the variable q than Eq. ffTUl) . Next we use the Feynman 
parametrization to evaluate this integral and we will see that the contribution to the d 
term is convergent. 
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We use the identity 



— — = 2 f dx I dy 



abc 



By setting 



o [a + (b — a)x + (c — b)y] 3 



a= (P -qf -m 2 K *, 

1, 2 2 

b = q - m K , 

c= (q — k) 2 — m 2 K 



and performing the change of variable 

q — q' + [(1 — x)P + yk] 

we have 

d 4 q' 1 



2 dx dy 
Jo Jo J 



(2vr) 4 (q' 2 + s + ie) 



4kM + [(l-x)P v + yk„]) 



+ 2( g ; t + [(1 - x)P, + yk,]){q' u + [(1 - x)P v + yk v \) 

(P-q' - [(1 - x)P + yk]) ■ (2k - 2P) + m 2 K , 



x 1 



rn 



K* 



with 



s = (1 — x)(xM 2 — m 2 K * — 2yP ■ k) — xm 2 K . 



(14) 



(15) 



(16) 



(17) 



However, in Eq. (1171) . all the terms that contribute to the d coefficient are finite. Hence, 
in the end, the evaluation of the amplitude of the type-b diagram is completely finite, and 
gives: 



1 m K *G v 
i=9biK* k& — 7^r 2P ■ ke A - e 



dx I dy 



V2P 
1 1 

32vr 2 s + ie 



-4(1 



x) 



Ay(l - x 



(xP - yk) -{k-P) 



m 



K* 



(19) 



In the derivation of Eq. f|T9|) from Eq. ffTTj) we have used the relation between the a and d 
coefficients given by Eq. ([6]). We have also used that [23] 



d A q' 



TX 



(q' 2 + s + ie) 



2 s 



is 



(20) 



and that terms with odd powers of q' vanish when performing the integration, by sym- 
metry reasons. It is also worth explaining a subtle cancellation which occurs between two 
logarithmic divergent pieces when deriving Eq. (|T9|) . as explained below: 
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In Eq. (fTTj) . apart from the terms which provide a finite contribution to the d-term, 
already considered in Eq. (fT9|) . there are two more terms which contribute to the d-term 
and which are logarithmically divergent. One of them goes as yk^q^ik — P) a . After 
the q' integration, this gives a term proportional to —yk^P u , since the q' v q' a gives a result 
proportional to g va . The other term goes as q' (1 — x)P u q f a (k — P) a and gives a term 
proportional to (l—x)k^P u after the q' integration, with the same proportionality coefficient 
as in the other case. However after doing the x and y integration these two terms give 
the same result but with opposite sign. Hence these two possible sources of divergent 
contribution to the k^Py cancel exactly among themselves. Therefore the expression of the 
amplitude in Eq. ( 1191) is totally finite. It is worth stressing again the power of the technique 
used here to evaluate the amplitude coming from the type-b loops since, despite starting 
from a loop quadratically divergent, we have been able to get rid of all the divergences in 
an exact way. 

At this point, it is worth noting that the numerical evaluation of the term proportional 
to 1/ m 2 ^* in Eq. (TT9"j) is about 5% of the other term. This term comes from the last factor of 
Eq. (fTTj) which essentially is due to the the p v p v /m v 2 part of the vector meson propagator. 
Hence, the 1 / my 2 terms can be safely ignored in the evaluation of the type-c diagram. This 
approximation is expected to be very accurate since, advancing some results, the type-c 
diagrams will be found to be very small compared to the type-b and hence the 1/my 2 is a 
small piece of a diagram contributing little to the radiative decay width. Nonetheless, we 
will include later on this uncertainty in the theoretical error analysis. 

Now we evaluate the amplitude corresponding to the type-c diagram, Fig. [lb) . 

We also need in these case the VV'j vertex that we get from gauging the charged vector 
meson kinetic term 

C = — F}. v F ixv (21) 

with F^ u = d^V u — d u V fl with the minimal coupling substitution <9 M — > + iqA^. 

After neglecting the 1/my 2 term of the vector meson propagator by the reasons ex- 
plained above, we have: 



1 m K *G v p 

_ c ~~ 1/2 K * ~sj2P~ €a€ 

with 

r d A q 1 1 1 

^ v J (2tt) 4 q 2 - m\, {q - kf - m 2 K , {P - q) 2 - m K 2 

x [2q u (2P - k - g) M - q ■ (2P - k - q)g„ v - (q - k)^(2P -k- q) v ] . (23) 

After keeping only the terms contributing to k^P^, doing the Feynman parametrization 
and using the relation of Eq. ([6]) the final expression of the amplitude coming from the 
type-c diagram is 



S 



1 

t c — J^gbiK'K^ 



2P 



keA ■ 6 



I dx dy 

J 



1 1 



[y(l -x)-3x + 2y + 1{24) 



V2P 



32tt 2 s' + is 



with 



s' = (l 



i)(iMj — m^; 2 — 2t/P • A;) — xm\ 



(25) 



Despite the smallness of the terms coming from the 1/my 2 term of the vector meson 
propagator in the type-b mechanism, it is worth mentioning a technicality regarding the 
cancellations of the divergences had we evaluated these terms in the type-c loop. If one 
keeps these 1/my 2 terms in the vector-meson propagators one obtains that the terms with 
1/my 4 do not contribute to the dk^P v term and there remains a logarithmic divergence 
proportional to 1/my 2 . This divergence should be expected to cancel had one included 
suitable tadpoles which could cancel the offshellness of the momentums of the vector-meson 
in the loops, in a similar way to what was shown in [24], where the factorization of the 
q 2 terms in the loop was justified. For the same reasons, this factorization was also used 
in [8]. Should one take this prescription here, the 1/my 2 terms would be also very small. 
In any case we will make a conservative estimate of the errors induced by considering these 
terms in one way or another. By knowing that the contribution to the width of the type-c 
loop diagram is one order of magnitude smaller than that of the type-b one, and that the 
changes found for the loop of type-b are of the order of 5%, an estimate of 10% error of 
the radiative width coming from these considerations is a safe estimate. 

Adding both type-b and type-c loops, we have: 



Another possible diagram with the photon directly emitted from the VPP vertex, which 
is needed to ensure gauge invariance of the set of diagrams, does not give contribution 
to the d coefficient since the vertices involved are both of the type e' ■ e, with e' either 
the vector or axial-vector meson polarization vector. Therefore there is no k momentum 
dependence either in the vertices or in the propagators and hence the integration cannot 
give contribution to the k^P u structure. 

Concerning the a\ — > 7r + 7 decay, the evaluation is totally analogous to the previous 
one. The amplitude from the K*K loops, when adding both type-b and -c mechanisms is: 




(26) 
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X f^-(l- X )[l + y( X P-yk).(k-P)/m 2 K .] 
\ S ~\ IS 

- -TT^W -x)~3x + 2y + 1}) (27) 

S ~\ IS J 

where one has to change by m ai in the evaluation of s and s' in Eqs. (1181) and (1251) . 
Note the relative minus sign in the terms with s' of Eq. (12?]) with respect to Eq. (1261) . This 
is due to the fact that, as we mentioned above, the b\ couples to the positive G— parity 
combination (K*K + K*K) while the a x couples to the negative G— parity combination 
(K*K-K*K). 

In the Gti case there is also the possibility of having p and n in the loops. This mechanism 
gives: 



t\Zl = 9a lP ,e^02P .ke A .e[dx[ dy^ 



y/2~p 71 Jo Jo y 327r 2 

x -x)[l + y{xP - yk) ■ (k - P)/m p 2 ] 

- 7^[2/(l -x)-3x + 2y+ 1}) (28) 
s + te J 

where one has to change by m ai , mx* by m p , and mx by in the evaluation of s 
and s' in Eqs. (TT8]) and (J2SJ) - 

The expression of the radiative decay width in the axial-vector meson rest frame is 

— > 

r < M -> = I2^ |T|2 < 29) 

where Ma stands for the mass of the decaying axial- vector meson and T is the sum of the 
amplitudes from the tree level and loop mechanisms removing the 6a • e factor. 

On the other hand, giving the large width of the axial-vector meson, particularly for 
the di, it is appropriate to fold the expression of the amplitude with the mass distribution 
of the axial- vector resonance. Hence the final amplitude is obtained from the expression 

i r (M A +2T A ) 2 dsA ( I ) r~r 

Ta-^p^M- 1 H— im\ .., ... ,. r(^I)e(^I-V4) (so) 

J(M A -2r A )z it \s A - M\ + iM A T A J v 

where G is the step function, Ma and V a are the nominal mass and total axial-vector meson 
width from the PDG [20] and s 1 ^ is the threshold for the dominant A decay channels. The 
errors quoted in the PDG for these magnitudes are taken into account in the error analysis. 
In Eq. ( 1301) . M is a normalization factor in the convolution integral obtained from the same 
integral as in Eq. (1301) setting T(- s J~sa) = 1. 

Once the formalism and the different vertices have been exposed, we are in a situation to 
address the possible contribution from mechanisms involving the n-ai mixing. The mixing 
of axial-vector and pseudoscalar mesons (or vectors and scalars) is possible through the 
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longitudinal component of the spin-1 propagator P^P V / P 2 [25-29]. In the Appendix we 
show that the diagrams corresponding to the present problem, involving this mixing, vanish 
in our formalism. 

3 Results 







T\ 

L t>1^7T7 


L a\ — >7T7 


nee level 




zo 






UJ 


16 






p 




650 




total 


76 


650 


loops 


type-b K*K 


30 


10 




type-b pir 




96 




type-b total 


30 


136 




type-c K*K 


0.14 


0.05 




type-c pit 




0.7 




type-c total 


0.14 


1.02 




K*K 


34 


8.7 




pit 




102 




total loops 


34 


137 


total this work 




210 ± 40 


460 ± 100 


experiment 




230 ± 60 [30] 


630 ± 246 [31] 



Table 1: Different contributions to the radiative decay widths. All the units are KeV. 

In table[T]we show the different contributions to the partial decay width coming from the 
different mechanisms considered in the calculation. The theoretical error in the final results 
have been obtained by doing a Monte-Carlo sampling of the parameters of the model within 
their uncertainty and considering the uncertainties in the couplings discussed above. Note, 
however, that we have no freedom in the theory once the relevant parameters (actually a 
subtraction constant) are fixed. To these errors we add in quadrature the 10% from the 
arguments used above concerning the 1/rriy 2 terms. 

From the results one can see that the tree level contribution for the a\ accounts for 
most of the decay width. However, for the b\ the tree level by itself only accounts for about 
1/3 of the experimental result, despite having two diagrams, <p and lo, that contribute to 
the tree level process. The smallness of the tree level contribution comes from the — v^/3 
and 1/3 factor of the and uj coupling to the photon in comparison to the factor 1 for the 
p case present in the a± decay and also to the fact that the a^pit coupling obtained with 
the chiral unitary model [8] is larger than the bi<j)ir and b\UJi\. 

Note also the constructive interference between the <fr and uj diagrams despite the coef- 
ficient of the V'y coupling having a different sign. This is so because the couplings of the 
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b\ to pre and ujtt have also relative different sign. These relative signs are also a genuine 
non-trivial prediction of the chiral unitary model of Ref . [8] . 

Regarding the loop contribution, the total loop results for the b\ and a\ decays have 
a comparable absolute value. In the 61 case it increases the decay rate to account very 
well for the observed experimental result, after interfering constructively with the tree level 
mechanism. Note that the most important contribution from the loops comes from the 
type-b mechanisms (see Fig. W?)- Particularly, for the a\ case the dominant contribution 
to the loops come from the pi\ loops. 

It is worth stressing the important role of the interferences among the different terms 
to give the final result. The interferences depend essentially on the sign of the couplings 
and the imaginary part of the loop functions. The values and relative signs of the AVP 
couplings are non-trivial predictions of the chiral unitary model of Ref. [8] and hence, the 
agreement of our calculated radiative decay widths with experiment gives support to the 
model of Ref. [8] and, hence, the dynamical nature of these axial-vector resonances. 

4 Conclusions 

We have studied here the radiative decays af — > 7r + 7 and bi — > 7r + 7 which had proved 
problematic before in several approaches. The novelty which allowed us to obtain a si- 
multaneous good description of both decays was the consideration of the ai and b\ axial 
vectors as dynamically generated resonances within the context of unitarized chiral pertur- 
bation theory. Because of that we found important loop contributions that were essential 
in reproducing the experimental values. Technically, it is particularly rewarding to see 
that, by invoking gauge invariance, the calculation is simplified and the relevant loops are 
shown to be convergent despite their large superficial degree of divergence. The nature 
of the resonances as quasibound states of meson and vector-meson states reverted into a 
loop contribution which provided a substantial contribution to the radiative amplitudes, 
particularly to the one of the 61 radiative decay. 

One might think that loop contributions could have been considered without resorting 
to the concept of a dynamically generated resonance, by simply taking the couplings of 
the resonance to their decay channels. However, for the important case of the b\ we found 
that the contribution of the ujtt loop vanished and the relevant contribution was coming 
from KK* and K*K which is a closed decay channel (up to the width of the states) 
and for which there is no valuable experimental information. The chiral unitary approach 
provides directly such couplings with definite signs since these states are a part of the 
building blocks of the resonances in the coupled channel approach. Similarly, with the use 
of a phenomenological Lagrangian, like the one of Ref. [13], one could get such couplings, 
but these are based on SU(3) symmetry which is actually broken when one generates 
dynamically resonances with a nonperturbative approach like the one in Ref. [8]. One 
example of relevance to the present case is that, with the phenomenological Lagrangian, 
the b\ — > 07r coupling is forbidden while in our case, the nonperturbative treatment of the 
problem, involving many iterative loops, generates a finite coupling that is dominant in 
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the tree level contribution of 61 — > 7r7 (see Tabled]). 

The fact that we obtain a good description of the two radiative decay rates for the first 
time provides support for the idea of the axial vector mesons as dynamically generated 
states within chiral dynamics. Other tests could follow as we get increased and more 
accurate information on the axial vector mesons, and the findings of the present work 
should serve to stimulate efforts in this direction. 



5 Appendix: Mechanisms related to the mixing of 
axial-vector and pseudoscalar mesons 

In addition to the terms discussed so far in this paper, we could have terms involving 
the mixing of the axial- vector and pseudoscalar mesons [25-29], through the longitudinal 
component of the axial-vector resonance. In our case this occurs with ai-w mixing, allowed 
by G-parity. The possible mechanisms involving this mixing in our scheme are given in 
Fig. [2j However, we shall demonstrate here that they vanish in our formalism. 



a) 



p 

b) 



Figure 2: Diagrams involving the a\-TC mixing. 



A free spin-1 massive meson propagator can be written as 



<r + 



pii pv 



<r + 



p/t pv 

-p3- 



+ 



pfipu 



P 2 -M 2 P 2 -M 2 ' M 2 P 2 ^ 

where in the second term of the equality a separation has been done in terms of a transverse 
part (—g^ u + P 2 ) and a longitudinal one (P^P"). Note that the pole of the particle 
appears only in the transverse part. 

In our formalism, the axial-vector resonance is dynamically generated from the VP 
interaction and is associated to the poles of the scattering matrix. In the Appendix B 
of Ref. [8] we made explicitly the separation into transverse and longitudinal part, with 
the result that the poles appeared only in the transverse part of the amplitude. There we 
found for Typ^v'P' 



Vb 



A1V 



pix pv 



Vc P^P V 
1-c P 2 



(32) 



where P is the total momentum of the VP system and e, e' , the polarizations of the two 
vector mesons. In Eq. fl32l) . c is very small compared to b and of opposite sign, such that 
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there are no poles in the longitudinal part. If we consider also ai-n mixing, we would have 
to add terms like in Fig. [3] to our VP amplitude. 




Figure 3: Extra contributions to the VP — > VP interaction involving the ai-ir mixing 

The loop function appearing in Fig.[3]has the structure J(P 2 )P^. The sum of terms in 
Fig. [3] renormalizes the longitudinal part of Eq. (1321 which is changed to 

Vc P»P V _ Vc(P 2 - m 2 ) P^P V 

1 _ c _ |gg [~P^~ ~ (l-c)(P 2 -m 2 ) -(3J 2 {P 2 ) P 2 (33) 

where m is the pion mass. The amplitude has the unphysical feature of providing a pole 
related to the pion pole (close to m 2 if (3 is small). The way to remove this unphysical 
behaviour of the longitudinal part is to demand that J(P 2 = m 2 ) = 0, which also appears 
in other formalism [28]. In other works [29] it is shown explicitly that the renormalized full 
vector meson propagator contains only one pole which does not show up in the longitudinal 
part. 

The contribution of the mechanisms of Fig. [2] in our formalism would have to be con- 
sidered through the a\ pole of the amplitudes of Fig. 0] 




Figure 4: Mechanisms of Fig. [2] in the dynamical formalism 



Diagram Fig. Hh) , with the photon emitted either from the pseudoscalar or the vector 
in the loop, is proportional to J{P 2 = m 2 ) and hence vanishes. 
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Diagram [4b ) is more subtle. The amplitude is proportional to 



£m(p) 



Vb 




) 



Vc P^P 



J{P 2 )P 



1-b 



+ 



/3J 2 (P 2 ) p2 
P 2 —m? 



1 - C 



e^P^P 2 ) 



Vc 



(34) 



1-c- 



f3J 2 (P 2 ) 



which has filtered the longitudinal part of the amplitude. Furthermore, the procedure we 
have followed to evaluate the coupling of a\ to rcy is equivalent to calculating the residue 
of the np — ► 7T7 amplitude and dividing by the a\ — > 7rp coupling [32]. The residue of 
Eq. ( 1341) at the resonance pole (6 = 1) is zero and hence the mechanism of Fig. 0b) also 
vanishes at the a\ pole. 
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